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Tunneling of a large spin via hyperfine interactions 
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We consider a large spin S in the magnetic field parallel to the uniaxial crystal field, interacting 
with N ^ 1 nuclear spins li via Hamiltonian H = -DS^ - H^S^, + AS-^^.^iIi with A < D, 
at temperature T. Tunneling splittings and the selection rules for the resonant values of are 
obtained perturbatively. The quantum coherence exists at T <C ASI while at T > ASI the coherence 
is destroyed and the relaxation of S is described by a stretched dependence which can be close to 
logt under certain conditions. Relevance to Mni2 acetate is discussed. 



PACS numbers: 75.45+j, 73.40Gk, 75.60.Lr 



The eflfect of nuclear spins on tunneling of a large 
spin S has beenj-ihe subject of intensive research in 
the last decadeJjTl The interest to this problem is two- 
fold. Firstly, nuclear spins are likely to significantly af- 
fect tunneling and coherence. Secondly, they represent a 
non-conventional dissipative environmentjlhat cannot be 
treated byjjhe Caldeira-Leggett method!] Models stud- 
ied to dateETH assume that tunneling of S is induced by 
those terms in the Hamiltonian, e.g., transverse field or 
transverse anisotropy, which doi-riat commute with, e.g., 
a longitudinal anisotropy term.QTj In that case nuclear 
spins interfere with tunneling by producing a bias field 
on the spin S and driving it off resonance. The rate 
of tunneling of S then decreases by a statistical factor 
that reflects the probability of the environmental spins 



arranging such that X^iLi = 0-^ 

In this paper we study a different situation when tun- 
neling is actually induced by the interaction of S with 
nuclear_apjiHS. This situation may be relevant to Mni2 
acetateliXllZl where the crystal-field terms responsible for 
tunneling are very small. We consider the following 
Hamiltonian 



H 



-DS^,~H,S,+ AS -Itot, Itot=^I», (1) 



where the electronic spin S ^ 1, Itot is the sum of iV ^ 1 
nuclear spins with magnitude /, and A <^ D is em effec- 
tive hyperfine constant. As an illustration we will keep 
in mind the example of Mni2 with 5 = 10, / = 5/2, 
TV = 12, D = 0.6 K, and A = 2 mK. We shall assume 
that the nuclear time Ti needed for 1^ to flip is larger than 
the tunneling time for S. The Hamiltonian conserves the 
magnitude /tot of Itot 7 thus dynamically it describes a 
two-spin problem within each subspace of fixed /tot- On 
the other hand, for an assembly of Mni2 molecules there 
is a distribution over the values of /tot, and all results 



should be averaged over this distribution. Contrary to 
previous approaches to problems with nuclear spins, we 
shall solve the above-formulated problem rigorously, us- 
ing the smallness of A in comparison with D. 

The longitudinal part of the hyperfine interaction, 
ASzItot.z splits the energy levels of the electronic spin 
which become 



= —Dm — HzTTL + Ammi, 



(2) 



where m and m/ are the projections of S and Itot on the 
z axis. The transverse part of the hyperfine interaction, 
(v4/2)(S'+/tot,- + S'-/tot,+), yields only small corrections 
to the energy levels in the wells in the second order of the 
perturbation theory which are of order /D and can be 
neglected. On the other hand, the transverse spin terms 
may cause transitions between degenerate energy levels 
on different sides of the barrier. 

In contrast to the problem of spin tunneling in the 
transverse field, here the projection of the total spin of 
the system S + Itot on the z axis is conserved. The value 
of the bias field for which resonant transitions occur, 
is defined by the conservation laws 



(3) 



and can be parametrized by the two quantum numbers 
k, kj 



H,{k,ki) = Dk 
where k and kj satisfy 



Akr 



m = — m — k, mj = ni + kj 
m'j = m — m' + mj = m + kj. 



(4) 



(5) 



Here, by definition, to < (see Fig. 1 of Ref. |l^ for 
an illustration of electronic resonances). Whereas the 



1 



number k is integer and labels the main electronic reso- 
nance, the number kj may be integer or half integer and 
describes its hyperfine splitting. For the zero-field reso- 
nance, k = kj — 0, the projection of the total spin of 
the system on the z axis is zero for both initial and final 
states. This only takes place, however, if both S and /tot 
are integer or half integer. If one of them is integer and 
another half integer, there is no resonance at zero field. 

For a fixed value of /tot, the possible values of kj are 
defined by the conditions —/tot < mi,m'j < /tot and 
satisfy 



k<ki <Iu 



(6) 



One can see that close to the top of the barrier, where 



and 



are close to each other, the number of values 



of kj is slightly less then 2/tot + 1- Since /tot < NI, the 
maximal number of these resonances is about 2NI + 1. 
On the contrary, for the unbiased ground-state resonance 
(to = —S, m' = S) there are only 2(/tot — 5) -I- 1 values 
of kj. For systems with NI < S, there are no resonant 
transitions between ground-state levels. 

Let us consider now the tunneling splitting Ae be- 
tween the resonant levels. Since the hyperfine constant 
A is much smaller than the anisotropy D, it can be 
calculaijed, jJiiih the help of the high-order perturbation 
theory.le'y N One has 

'^£m.mi;ni'.7n', = 2(to, TO/|//|to + 1, TO/ - 1) 

^ (to+1,TO/- l|//|TO-f2, TO/-2) 



^ rn + 1 , 771 / — 1 ^m/mj 

1 



-(to' — 1, TOj + 1|//|to', TOj). (7) 



Performing products in this formula, one arrives at the 
final result 



Ae, 



2{D + A) 



A 



I jS + toO!(^ - to)! (/tot - m;)!(/tot + mj)! 
(5-TO')!(5-l-TO)!(/tot+TO;)!(/tot-m/)!- ^ ' 

The part of this expression depending on the state of the 
nuclear subsystem, 



/nucl(/tot, — 



'(/t< 



fc7)!(/t, 



TO 



(/t< 



ki) 



, (9) 



has its maximum at kj in the middle of its interval, kj = 
— {m + m')/2 — k/2 [see Eq. (^], and the minimal value 
at the borders of this interval, kj = —/tot ~ ^ and kj — 
hot — to'. If the number of allowed values of kj is much 
greater than one, /nuci approaches a Gaussian (see Fig. 

For a given resonance (fc, kj) the splitting depends, 
apart from the initial and final electronic states, to and 
to', also on the value of /tot- Inverting Eq. (^), one ob- 
tains the range of possible values of /tot for a given reso- 
nance 




FIG. 1. Dependence of the level splitting Ae oc /nuci [see 
Eq. ^] on the nuclear quantum number ki in Eq. (^. It 
approaches a Gaussian for high spin values. 



max{TO' + ki, -TO - ki, (to' - to)/2} < /tot < NI. (10) 

The nuclear quantum number kj itself changes in the 
range ~NI ~ m < kj < NI — to', as follows from Eq. (||) 
if one replaces /tot by its maximal value NI. According 
to Eq. (||), one has m'j — —NI on the left border and 
TO/ = NI on the right border of this interval. For these 
border values of kj, there is obviously only one possible 
value /tot — NI satisfying Eq. ([lO|). On the other hand, 
for kj — k/2 the number of values of /tot attains its 
maximal value of about NI — to' — k/2. In addition, each 
value of /tot can be built up of individual nuclear spins 
in a number of different ways, if /tot < NI. The latter 
shows that the hyperfine resonances near the center of 
the splitted electronic resonance, kj = k/2, are much 
stronger than those for fc/ away from the center. It is, 
however, difficult to parametrize the strength of these 
hyperfine lines since the resonance dynamics is rather 
complicated, as will be seen shortly. 

The dependence of the level splitting on /tot for a given 
hyperfine resonance results in the decoherence of tunnel- 
ing. If electronic spins are prepared in the state to, then 
only molecules with to/ satisfying the resonance condi- 
tion of Eq. (H) will take part in the process. Among those 
with a given to/, there are molecules with different pos- 
sible values of /tot, which will show oscillations between 
the degenerate states with different frequencies. If in the 
initial state nuclear spins are in equilibrium, the proba- 
bility of finding the electronic spin in the initial state to 
depends on time according to 



Pmit) = 1 



^AA(/tot)[l-cos(A£t)], (11) 



where Zj = sinh{[l + 1/(2/)]^}/ sinh[C/(2/)] with ^ = 
Ami /T is the partition function of an individual nuclear 
spin, /tot satisfies Eq. (|l^), Ae = Ae,n,m,^r,i' ,m'j{Itot) is 
given by Eq. (|^), and M{Itat) is the number of ways to 
combine /tot- For temperatures larger than the charac- 
teristic nuclear temperature, ^ <C 1, the statistical factor 



2 



in front of the sum simplifies to (2/ + Under typ- 

ical conditions, only a small part of the molecules have 
a given value of mj and thus take part in the resonant 
tunneling, thus the dynamical term in Eq. ( pi] ) is much 
smaller than one. This is a mechanism of tunneling re- 
duction due to nuclear spins pointed out in Refs. ^ and 
^. The rest of molecules which do not have the required 
value of mi for a given resonance are "frozen in" , and 
they should wait for a longer time Ti required for the 
nuclear spins to relax. The study of the latter is beyond 
the scope of this paper. 

The number of realizations Af{Itot) in Eq. ( pT| ) can 
be computed recurrently. If the total spin of a system 
of N nuclei is /tot, the total spin of its subsystem of 
N — 1 nuclei L' * assumes the values IL, 



< 



min{/tot + I, (N — 1)1}. Thus for the number of realiza- 
tions Af{Itot, N) one can write. 



min{/tot+/,(JV-l)/} 



/tot 



MilU,N~l). (12) 



The initial condition for this recurrence relation is 
A/'(/t'ot,2) = 1 for < < 21. The quantity A/'(/tot) 
obeys the normalization condition 



NI 

E 

/tot=frac(Ar/) 



(2/tot + l)AA(/tot) - (2/+1) 



N 



(13) 



For NI > 1, the quantity (2/tot + l)A/'(/tot)/(2/-H 1)^ is 
the high-temperature distribution function of the magni- 
tude of /tot and it is well approximated by Airl^^^F^ltot), 
where /"(Itot) is a normalized Gaussian function with re- 
spect to the three components of Itot -t^ Thus one has the 
asymptotic form for NI ^ 1 



P(/t< 



27r/t, 



AA(/tot) ^ 

(2/+ 1)^ ~ (27rcr/)3/2 



exp 



2<TI 



(14) 



where c/ = (7V/3)/(/ -I- 1). It has a maximum at 
/tot = ^ which is about 6 for Mni2 (/ = 5/2, N = 12). 
Fig. ^ shows an agreement between the exactly computed 
P(/tot) and its Gaussian approximation for / = 5/2 and 
iV = 12. This agreement improves for higher values of 
NI. 

Due to the summation over /tot in Eq. ( ITTI) and the 
dependence of the splitting Ae on /tot, the coherence of 
tunneling is destroyed. On the other hand, tunneling 
oscillations are damped due to the finite linewidth of the 
electronic levels 7™. In the overdamped case 7™ ^ Ae 
the relaxation of the initially prepared state m to the 
matching state m' for ^ = Ami /T ^ 1 is described by a 
sum of exponentials (cf. Ref. [I9| ) 

p^[t) = l-^p(/tot)i[l-exp(-r„,„,,fc,,,,„,0], (15) 



where /tot satisfies Eq. (nC 
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FIG. 2. The normalized number of realizations of the 
total nuclear spin hot, P{hot) = Af{Itot)/i'2I+ 1)^: exact 
and Gaussian, see Eq. (pl|). 
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FIG. 3. Stretched relaxation curve, Eq. (^_8|), for the tun- 
neling resonances due to nuclear spins. The sets of parame- 
ters labeled by "The same" are variations of the basic set of 
parameters in the upper right corner. 



(Ae fin 



1mm' + (^1 



(16) 



7mm' = 7m + 7m', and e,n,mi " £m' ,m'^ = (m - m'){II - 
Hk^ki). As we have pointed out above, the process de- 
scribed by Eq. ( p^ does not lead to the full relaxation 
since only a small fraction of the systems possesses the 
required value of the projection mj = m' + ki of the total 
nuclear spin. The probability pm decreases due to this 
process by 

oo 

^P™ = ^ E -I J ^(^tot)rf/tot 



1 1 



2 y'27r(T/ 



exp 



tot.min 
2(7/ 



(17) 



where Eq. (^ was used and /tot,min is defined by Eq. 
(|To|). In particular, for the central hyperfine line of an 
electronic resonance, kj = fc/2, one has /tot.min — {fn' — 



3 



m)/2. For resonances between low- lying electronic levels 
this value is big and hence Apm is small. An extreme 
case is the resonance between m = —S and rn' = S for 
which — 0.008 for the Mni2 set of parameters. 

One should note that the relaxation described by Eq. 
( p^ ) is not a single exponential but a sum of different 
exponentials with a faster rate and smaller amplitude 
(larger /tot) and slower rate and larger amplitude (smaller 
/tot)- To illustrate this fact we have plotted in Fig. ^ the 
relaxation function 



R{i) = 



Ez,„,^^atot)exp[- 



(18) 



where /nuci is given by Eq. (||) and i is the scaled time 
including all factors which do not depend on /tot and kj. 
One can see that for the chosen values of the parameters, 
the relaxation function R{t) is stretched for about 14 
decades, whereas a single exponential practically decays 
within two decades in time. To analytically clarify the 
large-time behavior of R(t) one can apply the Stirling 
formula to Eq. (||) to obtain 



f2 ^ (pL \2{rn'-m) _._ 



1 



f m'+ki \ 



/ m+ki \ ' 

\ hat J 



m-i-ki 



(19) 



Here the second factor is only weakly dependent on /tot- 
For large m' — m, the dominant dependence on /tot is 
give n by the first factor. Substituting Eq. ( p^ ) into Eq. 

(|l8|), replacing the sum by an integral and introducing 

1 

/tot = /tot f^t™'-""' as a new integration variable, one 
obtains from Eqs. (Ill) and (p^ the large-time behavior 



R(i) cx r 



(20) 



Since m! — m can be large, e.g., about ten, the exponent 
can be rather small. Therefore Eq. (^^ can be approxi- 
mated by : 



R[i) cx 1 



■Inf, 



(21) 



which shows that the relaxation, although being a power 
law, may be well approximated by a logarithmic time 
dependence. We remind the reader that m' — m is always 
positive. 

For temperatures much lower than the nuclear tem- 
perature, i.e., for ^ = Ami /T ^ 1, the coherence is 
restored, since the nuclear spins fall in the lowest-energy 
state mi = /tot = NI^ and there is only one term in Eq. 
(pi]). This term corresponds to m = —S, m' = S — k, 
mi = NI, and m'j ^ NI -2S + k, i.e., ki = NI - S + k 
(see Fig. ^). Consequently, the low-temperature reso- 
nances are characterized by only one number k: 



H,{k) = A{NI -S) + iD + A)k. 



(22) 



In the case of Mni2 the corresponding quantum relax- 
ation due to nuclear spins can become observable at 



NI 



1 <^ 2 




NI-2S 



FIG. 4. Tunneling between electronic ground states ±S 
at zero temperature. The z projection of the total spin is 
conserved. To conserve energy, the field Hz — Aki with 
ki = NI — S should be applied [see Eq. (|)]. 

TABLE I. Tunneling splittings due to the nuclear spins for 



the zero-temperature resonances m — 


-S,m' = S -k for the 


Mni2 set of parameters. 




k 


Ae, K 


10 


2.273148E-18 


11 


4.679795E-16 


12 


8.019877E-14 


13 


1.126574E-11 


14 


1.268525E-09 


15 


1.107789E-07 


16 


7.129369E-06 


17 


3.102200E-04 


18 


7.706470E-03 


19 


6.928205E-02 



k > 13, as is illustrated by Table 1 that lists tunnel- 
ing splittings for 10 < k < 19. Note that for large 
values of fc, tunneling due to hyperfine interactions in 
Mni2 should dominate over the effect of other small non- 
diagonal terms unaccounted for in Eq. (|^). 

In fact, the simplest model used in this paper cannot 
be quantitatively accurate for Mni2 since each molecule 
of Mni2 contains eight Mn atoms with the spin value 2 
and four Mn atoms with the spin 3/2. The latter couple 
ferromagnetically within each group, and the two groups 
couple ferrimagnetically with each other to build the to- 
tal spin iS* = 10 (see the details in the recent Ref. |2^). The 
hyperfine Hamiltonian of a Mni2 molecule contains two. 
hyperfine constants A and A' for each type of Mn atoms.Q 
Thus nuclear spins are splitted into two groups which be- 
have dynamically as two effective "giant" nuclear spins, if 
the nuclear relaxation is neglected. The resulting three- 
spin model (one electronic and two nuclear spins) is more 
complicated than the simplified two-spin model consid- 
ered above, although tunneling splittings between differ- 
ent resonant states can be obtained by the same method, 
cf. Eq. (^. The resonant values of the bias field Hz are 
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parametrized by four quantum numbers instead of two M. I. Katsnelson, V. V. Dobrovitski, and B. N. Harmon 

in Eq. and hence there are much more hyperfine hues Phys. Rev. B 59, 6919 (1999) 

for each electronic resonance. We do not try to work 
out this more reahstic model here since the calculations 
should be rather cumbersome and the accurate values 
of A and A' are unknown. This, however, does not af- 
fect our qualitative conlusions derived from a simplified 
model. 

In conclusion, we have studied tunneling and relax- 
ation of a large spin, induced by the hyperfine interac- 
tions. This model can be relevant to Mni2 acetate in the 
absence of the transverse field. Our main findings are 
these. There must be two distinct temperature regimes, 
above and below the hyperfine temperature, T^f = ASI, 
which for Mni2 is about 50 mK. At T ^ Thf the mag- 
netic relaxation is stretched and close to logarithmic due 
to the contribution of a large number of nuclear states. 
At T <C T/i/ only one nuclear state contributes to the 
magnetic relaxation at each resonant field, and the co- 
herence of tunneling is restored. 
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